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ABSTRACT 

To constrain the giant pulse (GP) emission mechanism and test the model of Lyutikov (2007) for 
GP emission, we have carried out a campaign of simultaneous observations of the Crab pulsar at 7-ray 
{Fermi) and radio (Green Bank Telescope) wavelengths. Over 10 hours of simultaneous observations 
we obtained a sample of 2.1x10^ giant pulses, observed at a radio frequency of 9 GHz, and 77 Fermi 
photons, with energies between 100 MeV and 5 GeV. The majority of CPs came from the interpulse 
(IP) phase window. We found no change in the GP generation rate within 10— 120 s windows at lags 
of up to ±40min of observed 7-ray photons. The 95% upper limit for a 7-ray flux enhancement in 
pulsed emission phase window around all GPs is 4 times the average pulsed 7-ray flux from the Crab. 
For the subset of IP GPs, the enhancement upper limit, within the IP emission window, is 12 times 
the average pulsed 7-ray flux. These results suggest that GPs, at least high-frequency IP GPs, are 
due to changes in coherence of radio emission rather than an overall increase in the magnetospheric 
particle density. 

Subject headings: Crab pulsar. Giant Pulses, Fermi 



1. INTRODUCTION 

The Crab pulsa r was discovered by 
iStaelin fc Reifenstein IIJ in 119681 by its remarkably 
bright giant pulses (GPs). Giant pulses are short (from 
few ns t o few ^ts), sporadic bursts of pulsar radi o 
emission (jPopov fc Stappersl[2007l : iHankins et an[200l . 
The nature of GPs is far from being clear and even 
the precise definitio n of g iant pulse had not yet been 
given (Knight ct aQ [20M ). GPs generally occur only 
in certain narrow ranges of pulse phase that are often 
coinc ident with puls e s seen at X-ray an d 7-ray ener- 
gies (jLundgrenl 119941 ). iPopov et all (|2006[ ) propose that 
all radio emission from the Crab (except for that in the 
precursor) is composed entirely of GPs, consistent with 
the alignment of the GP and high-en ergy components 
seen in other pulsars ex hibiting GPs (jCusumano et ahl 
[2OOI Ifcfight et al.l[2006h . 

The Crab pulsar shows pulsed emission across the 
entire electromagnetic spectrum (see Fig. [1] left), re- 
flecting different radiation processes in the pulsar mag- 
netosphere — from coherent curvature or synchrotron 
(radio) to incoherent synchrotron (optical and X-ray) 
and incoherent curvature (7-ray) radiation. Simi- 
lar to other sporadic variability phenomena seen in 
pulsar radio emission, represented by nulling pul- 
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sars f e.g. iHerfindal fc RankinI l2009f l. intermittent pul- 
sars (iKramer et al.~l2006[l. an d rotating radio tran- 
sients iMcLauSIinieraiJlool), GP emission could be 
due to changes in the coherence of the radio emission, 
variations in the pair creation rate in the magnetosphere, 
or changes in the beaming direction. If the GP phe- 
nomenon is due to changes in the coherence of the radio 
emission mechanism, then one would expect little cor- 
relation of the radio GPs with the high-energy emission. 
However, if the GPs are due to changes in the actual rate 
of pair creation in the pulsar magnetosphere, one would 
expect an increased flux at high energies at the time of 
the GPs. Similarly, because the radio GP and 7-ray com- 
ponents are aligned, one expects that they come from the 
same place in the pulsar magnetosphere. Therefore, if a 
GP occurs from a beam direction alteration, one would 
e xpect to also see an inc rease in the high-energy flux. 

iLundgren et al.l ()1995f ) previously attempted to carry 
out simultaneous radio/7-ray observations (50 — 220keV, 
the energy range of CGRO/OSSE) and correlate times 
of arrival of GPs at 800 and 1300 MHz with 7-ray 
photons. Their upper limit on the 7-ray flux in- 
crease concurrent with radio GPs was < 2.5. Later, 
iRamanamurthv fc ThompsonI ()1998f ) correlated the same 
set of GPs with EGRET photons of energy > 50 MeV, 
placing an upper limit on concurrent 7-ray flux of 4.6 
times the average Crab flux. This suggested that the 
GP mechanism is largely based on changes in coher- 
ence and not changes in pair p roduction rates or beam- 
ing. Yet, iShearer et al.l ()2003[ ) performed simultaneous 
radio/optical observations of the Crab pulsar and found 
a weak correlation, i.e. that optical pulses coincident 
with radio GPs were on average 3% brighter than others. 
This observation suggested that the GP emission mech- 
anism, whatever its nature, includes small variations in 
magnetospheric p article density. 

iLyutikovl ()2007l ) proposed a more specific, quantitative 
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model of GP emission in which Crab GPs are generated 
on closed magnetic field lines near the light cylinder via 
anomalous cyclotron resonance on the ordinary mode. 
During emission of a photon, an electron undergoes tran- 
sition up in Landau levels. The ene rgy is supplied by 
the parallel motion (jGinzburgI I1985D . The application 
of anomalous cyclotron reso nance to pulsar radio emis- 
sion has been discussed by iLvutikov et al.l (|1999[ ) and 
IMachabeh &:"Ui5^ ([l97l . 

One clear prediction of this model is that radio GPs 
(at least those at radio frequencies > 4 GHz) should be 
accompanied by 7-ray photons, as the high energy beam 
is expected to produce curvature radiation at energies 
^ h^^Vl ^ 0.1-100 GcV, depending on the exact value of 
the Lorentz factor 7. These energies fall into the energy 
range of the Fermi mission, and so this hypothesis can 
also be tested through high-frequency radio observations 
concurrent with Fermi. 

The Fermi Large Area Telescope (LAT) , with its large 
effective area, broad field of view, and superior angular 
resolution is a perfect tool for testing the Lyutikov theory 
and investigating the possible correlation between GPs 
and 7-ray photons in general. For the radio observations, 
using the 100-m Green Bank Telescope (GBT) allows one 
to record a very large number of GPs within a reasonable 
observing time, even at frequencies above 4 GHz. Thus, 
a thorough study of the correlation between high-energy 
7-ray photons and high-frequency GPs is possible. 

In this paper, wc present the results of simultaneous 
GBT/LAT observations of the Crab pulsar. To probe 
the level of correlation between GPs and 7-ray photons, 
we used two main approaches. First, we searched for 
a "burst correlation" by examining whether GPs cluster 
near 7-ray photons in time. Second, we analyzed whether 
the average 7-ray flux of the pulsar increases within the 
pulse phase windows where single GPs are produced. 

In Sections [2] and [3] below we describe the radio obser- 
vations and Fermi data used in this analysis. Section 2] 
discusses the influence of the interstellar medium on the 
observed GP sample. Wc describe the correlation analy- 
sis between radio GPs and Fermi photons in Section [5l 
and conclude in Section [6l 

2. RADIO OBSERVATIONS 

The radio observations were carried out during 12 ob- 
serving sessions in September-October, 2009 with the 
GBT, using the new Green Bank Ultimate Pulsar Pro- 
cessor Instrument (GUPPI) at a central frequency of 
8.9 GHz, in incoherent dedispersion mode. The total 
bandwidth of 800 MHz was split into 256 frequency chan- 
nels, and the total intensity was recorded with a sampling 
interval 2.56— 3.84 /is. Total observing time was ~ 26hrs 
or ~ 3 X 10^ pulsar periods. 

The raw data from every session were dedispcrsed with 
the current DM of the Crab pulsar^ using PRESTO pack- 
age^, and searched for all single-pulse events with signal- 
to-noise ratio (S/N) > 7. Since GPs from the Crab pulsar 
do n ot have any established lower limit on peak flux den- 
sity (jPopov et a l. 2006), we picked up initial threshold of 

* The DM was 56.8005 pc cm"'' for September and 
56.8109 pc cm~^ for October, from the Jodrell Bank Crab pulsar 
mo nthly ephemeris: http : //www . jb . man . ac . uk/pulsar/crab . html , 

® |http : //www ■ cv . nrao ■ edu/$\sim$sraiisom/presto/ 1 
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Fig. 1. — Left. Avera ge profile of the Crab pulsa r from radio to 7- 
rays, from the paper of lMoltett fc Mankind l l 19961 ). Right. Average 
Crab pulsar radio profile for one out of two sub-sessions with the 
GBT on Sep 25, 2009. 



S/N > 7 in order not to contaminate our sample by nu- 
merous spurious detections on noise. Each event was as- 
signed a width, found by averaging the dedispersed time 
series with different numbers of samples and finding the 
number that resulted in a peak in S / N. The list of event 
times was put into TEMP02 format (jHobbs et al.ll2006[) 
and converted to the barycentric reference frame for the 
correlation analysis with Fermi data. Times of arrival 
(TOAs) were corrected for delay due to propagation in 
the ionized interstellar medium (ISM). 

Estimated timing errors due to an inaccurate DM are 
less than our time resolution, assuming that DM varies 
smoothly and that between observing sessions the change 
in DM is less than the change over two months. For 
DMoct — DMscp = 0.0104pccm"'^, timing errors are 
about 0.5 /is. 

Fig. [1] (right) shows the average pulse profile (top) of 
the Crab pulsar at 8.9 GHz together with the subintegra- 
tions from one of the two sub-sessions on Sep 25, 2009, 
which had the highest rate of GP detection of all 12 ses- 
sions. The interpulse (IP) and high-frequency compo- 
nents (HFCs) are clearly seen, with the weak peak after 
HFC2 being the main pulse (MP). However, on Sep 25 
the pulsar was the brightest, and during other sessions 
the average profile was less prominent. During some ses- 
sions we did not accumulate a detectable average profile 
at all. 

The system equivalent flux density (SEED) is mostly 
determined by the Crab Nebula. Flux densities for the 
Crab Nebula were calculated with the relation S{f) = 
955 X (//GHz)-0-27 Jy (jCordes et all [Mil , accounting 
for the fact that at 8.9 GHz the sohd angle of the GBT 
beam covers only 25% of the area occupied by the nebula. 
Wc estimate a SEED 1.3/^At/l/xs Jy, or about 0.7 Jy 
for our most common sampling time, or At, of 3.2 /iS. 
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TABLE 1 

Summary of observational parameters and GP/7-RAY outcome for each observing 
DATE. Columns include (from left to right): date of observation, time resolution, 
system equivalent flux density (SEFD), total duration of radio observations and 
the time simultaneous with Fermi, number of giant pulses, Nqps: detected during 

THE whole observing SESSION AND DURING THE TIME SIMULTANEOUS WITH Fermi, NUMBER 

OF 7-RAY PHOTONS, . 



Date 


At 


SEFD 


Radio time 




A^GPs 










Total 


With Fermi 


All 


With 




(2009) 


(/is) 


(Jy) 


(min) 


(min) 




Fermi 




Sep 12 


2.56 


0.81 


87.1 


28.9 


139 


8 


5 


Sep 14 


3.20 


0.73 


165.7 


62.1 


4375 


1834 


10 


Sep 16 


3.84 


0.66 


99.3 


30.6 


98 


27 


4 


Sep 19 


3.20 


0.73 


118.5 


54.3 


6957 


1830 


5 


Sep 20 


3.20 


0.73 


110.5 


32.2 


1846 


384 


2 


Sep 21 


3.84 


0.66 


55.1 


31.5 


27 


16 


2 


Sep 22 


3.20 


0.73 


147.7 


68.9 


1256 


603 


5 


Sep 23 


3.20 


0.73 


164.5 


82.0 


10520 


5078 


10 


Sep 24 


3.84 


0.66 


55.3 


22.5 


38 


37 


1 


Sep 25 


3.20 


0.73 


236.5 


130.8 


14320 


10014 


13 


Sep 28 


3.20 


0.73 


72.3 


48.1 


34 


9 


5 


Oct 25 


3.20 


0.73 


157.6 


41.3 


3164 


1261 


15 


Total 






1470.0 


633.1 


42774 


21092 


77 



Preliminary analysis of all events with S/N > 7 re- 
vealed that GPs appear mostly in the MP and IP phase 
windows. Only for 2 observing sessions, namely on Sep 
25 and 28, were there several GPs detected in the HFCs, 
all w ith S/N close to 9. T his is to some extent surprising, 
since lJessner et aD ()2005[ ) observed the Crab pulsar with 
similar parameters and timespan, and found about 120 
GPs in HFCs versus 350 GPs in the MP and IP. How- 
ever, their threshold peak flux density of 25 Jy was much 
higher than our threshold of about 6 Jy, suggesting that 
the GPs from HFCs are rarer but brighter. 

Some of our sessions were heavily contaminated with 
broad radio frequency interference (RFI) pulses with typ- 
ical S/N < 10. Therefore, we analyzed only events with a 
peak flux density exceeding 8.1 Jy (S/N = 10 for the ses- 
sion with smallest sampling time) and which arrived in 
the MP or IP phase windows. Additionally, we excluded 
all events with width larger than 30 sample intervals, as 
being presumably caused by RFI. These cuts resulted in 
the selection of more than 40000 GPs. Comparing the 
number of single pulses above 8.1 Jy and narrower than 
30 samples in and out of the pulsed emission phase win- 
dows, we can estimate the proportion oi false GPs in our 
final data set to be less than 0.001%. 

The summary of observations is given in Table [T] For 
each observing date the listed columns are: time resolu- 
tion, SEFD, total duration of radio observations and the 
time simultaneous with Fermi, number of giant pulses, 
detected during the whole observing session and during 
the time simultaneous with Fermi, and the number of 7- 
ray photons selected for further analysis (see Section [3|). 

3. FERMI DATA 

We extracted "Diffuse" class events with energies from 
100 MeV to 300 GeV from the Fermi database, concur- 
rent with each radio observation. Photons with zenith 
angles greater than 105° were excluded to eliminate the 
7-rays generated in the Earth's atmosphere. 

We selected only photons in Good Time Inter- 
vals (GTIs) within an angle 9 < Max(6.68 - 
1.761og(£:/1000 MeV), 1.3)° of the radio pulsar position 



(|Abdo et al.ll20ldl ). Photon arrival times were converted 
to the Solar System barycenter and assigned phases with 
the TEMP02 f ermi plugin. The timing accuracy of the 
Fermi LAT is better than 1 ^s (.Abdo et al., .20101 ) . LAT 
dead time per event is less than lOOyits^^, which is less 
than 3% of pulsar rotational phase. Over the course of all 
radio observations we accumulated 10.5 hours of Fermi 
data within GTIs, resulting in 77 photons with energies 
above 100 MeV (see Table ft]). 

Fig.lHgives a quick visual summary of our simultaneous 
observations, showing Fermi photons and radio GPs ver- 
sus observing time for each session. The distribution of 
number and energy/peak intensity of photons/CPs with 
respect to pulsar rotational phase is shown in Fig [3l As 
reported earlier, 7-ray and radio emission windows are 
aligned. 

4. PROPAGATION EFFECTS FOR GIANT PULSES 

At high frequencies . GPs are strongly a ffected by inter- 
stellar scintillations (jCordes et al.l[20M ). which change 
their apparent rate and peak intensities. Careful treat- 
ment of ISM effects is crucial for a proper correlation 
analysis. Unfortunately, we did not make direct measure- 
ments of typical ISM diagnostic parameters, such as scin- 
tillation bandwidth and pulse broadening time. Instead, 
we use the scintillation timescales scaled from other fre- 
quencies. These values give only a rough estimate of 
scintillation timescales at 8.96 GHz, since for the Crab 
pulsar the main contribution to scintillations is made by 
the turbulent and quickly changing Crab Nebula. Scin- 
tillation parameters for the Cra b have been shown to be 
strongly variable with time (see lCordes et al.l ()2004D and 
references therein). 

The refractive interstellar scintillation (RISS) 
timescale, triss, can be scaled using the j/"^ -^ de- 
endence derived fro m a five frequency data set 
Rickett fc Lvnelll990l ). At 8.96 GHz, triss is about 
80 minutes, thus roughly matching the observed day-to- 

The dead time was taken from Fermi Technical Handbook, 
http : / /f ermi . gsf c .nasa. gov/ ssc/proposals/manual/ 
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Fig. 2. — Time series of radio GPs and Fermi photons during 12 observing sessions. X-axis — time from the beginning of each session, 
in minutes. Y-axis (left) - peak flux density of radio GPs. Y-axis (right) - energy of 7-ray photons. Both scales are the same for each 
observing session. The yellow shaded regions mark the time when we actually were recording radio data and the blue are the Fermi Good 
Time Intervals. The observing date is given in the upper right-hand corner of each subplot. For the observing session on Oct 25 two 
photons came within short time interval, so their markers overlap and one can see only 14 photons, instead of 15. 
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Fig. 3. — Left: histograms of GPs for all radio observing time (bottom) and Fermi photons during the simultaneous time (top). For 
illustrative purposes, scaled radio and 7-ray profiles (grey) are shown with arbitrary offset along y-axis. The scaled radio profile is from the 
second scan of GBT session on Sep 25, 2009 (as shown in Fig.[T]|, and the 7-ray profile is the Fermi profile accumulated during Sep-Oct, 
2009. Right: distribution of peak fiux density of GPs and energy of 7-ray photons over pulsar rotational phase. Also, there is shown the 
scaled radio profile as for histogram on the left. 
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Fig. 4. — Left: cumulative peak flux distribution of GPs (MP and IP together) for different observing sessions and 2 subsessions of 



Sep 25 with apparently different GP rates. Poisson errors are also shown. Right : 
observing session. 

day variation of the GP rate (see Fig. [2]). Also, for the 
longest, 4-h session, note the change of GP rate between 
two 2-h scans in Fig. [2] (observing settings were the 
same for bot h scans). 

Following iCordes et al.l (|2004[ ). for calculating the 
diffractive interstellar scintillation timescale, tdisSj we 
adopted the thin screen model with a Kolmogorov spec- 
trum of irregularities and reference pulse broadening 
time Td = 0.5 ms at 0.3 GHz. At 8.9 GHz, with band- 
width of 800 MHz, this gives us a scintillation strength 
u = 8.9, well into the strong scintillaton regime, and 
T"Diss = triss/m^ - 9min. 

The variation of the GP rate within each observing 
subsession was estimated by autocorrelating the rate of 
GP emission in 10-second bins. The autocorrelation 
analysis shows two shorter GP rate variability time scales 
of ~20 minutes and 2—4 minutes. Both of them agree 
fairly well with the DISS estimations, considering all the 
uncertainty in the scintillation parameters. On the other 
hand, there is no evidence against the hypothesis that at 
least one of these timescales is due to intrinsic GP rate 
variability. 

The following analysis assumes that observed day-to- 
day variation ( or, in case of Sep 25, variation between 2 
subsessions) of the GP rate and mean intensity is caused 
by RISS. If intrinsic GP rate and mean intensity are 
constant on timescales larger than 90 min, then it is rel- 
atively easy to make a GP sample corrected for refrac- 
tive scintillation. We accomplish this by estimating the 
amount of RISS intensity damping on each observation 
session with respect to the session with the highest GP 
rate. Then, we multiplied the intensities of all pulses in 
each separate session by those amounts, and threw out 
all GPs below a threshold, common for the corrected GPs 
over all sessions. 

A simple way to calculate the intensity variation due 
to RISS would be by comparing mean profiles of pulsed 
emission accumulated during each session. However, at 
these frequencies, our observations were not sensitive 
enough to accumulate the normal Crab pulse profile ex- 
cept on one or two sessions where scintillations caused 
a boosting of the average flux density of the pulsar. In- 



X values versus the RISS damping coefficient, k, for each 



TABLE 2 

RISS CORRECTION COEFFICIENTS k FOR EACH OBSERVING SESSION. 
Session RISS correction coefficient k 



Sep 12 


6.4 


3.1 


Sep 14 


2.2 


21.7 


Sep 16 


7.2 


4.7 


Sep 19 


1.4 


3.8 


Sep 20 


2.7 


3.9 


Sep 21 


8.5 


3.3 


Sep 22 


3.4 


6.6 


Sep 23 


1.3 


4.8 


Sep 24 


7.7 


2.5 


Sep 25/1 


5.0 


2.0 


Sep 25/2 


1.0 


0.0 


Sep 28 


8.5 


1.7 


Oct 25 


2.5 


5.2 



stead, we compared the intensity distributions of GPs 
between sessions. If the change in rate and mean inten- 
sity of GPs on timescales of a few hours is due to RISS, 
then the peak intensity distributions for each day should 
have the same shape, but with different values of peak 
flux density. The distributions in Fig. 2] (left) show that 
this assumption is basically correct. 

As the reference session, we picked the one with the 
highest rate of GPs, the second subsession of Sep 25 
(from now on called "0925/2" or "reference session"). 
For each day and for the two subsessions on Sep 25 sep- 
arately, we determined the RISS damping coefficient, k, 
by minimizing x^- 



xl 



1 



[iV„f (/,) - N{I,/k)]' 



Nuns - 1 Y 



(1) 



'N{U/k) 



where N{Ii) is the number of GPs with peak intensity 
higher than Ii, per hour of observation. A^bins is the 
number of bins in the distributions being compared, and 

<JN{ii/k) = \/ N{Ii/k)/Tiirs is the Poisson error in each 
bin if Tiii-s is the duration of the session in hours and 
assuming that the energy of each GP does not depend 
on the energy of the preceding one. 
The Xk curves are plotted in Fig. |4] (right). All but one 
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Fig. 5. — The distribution of time lags between GPs and 7-ray photons for 2-min bins, for all GPs > 8.1 Jy (21000 GPs; left) and the 
RISS-corrected GPs (180 GPs; right). Real Fermi photons (blue line) are contrasted to the mean and 95% percentile on the pool of 
simulated data sets (grey and black). The fact that the distribution for the real data lies within 95% of the simulated ones indicates no 
apparent change in GP generation rate on 2-min timescale with any possible time lag up to ±40 minutes with respect to the 7-ray photons. 
The maximum lag value, 40 min, corresponds to the size of the largest GTI window. All other bin widths (down to 10 s) give the same 
result. 



have a sharp minimum oi xi ^ 10) indicating reasonable 
fits, at fcs between 1 and 10. The only outlier is the 
session of Sep 14th, with Xmin ^ 21.7, which has an 
abnormal excess of intrinsically strong GPs (see Fig. [2]). 
These GPs do not exhibit any other peculiar properties, 
other than relatively high peak flux density. 

Both RISS coefhcients and the corresponding Xmin 
for each day are listed in Table [2j We corrected the 
session GPs by multiplying their flux densities by k 
and set the intensity threshold for corrected pulses as 
8.1 Jyx max(fc) — 69 Jy (8.1 Jy was our initial threshold, 
see Section [2l). 

Thus, we effectively selected only those pulses which 
would have had peak flux density larger than 8.1 Jy 
if they were observed during the session with highest 
RISS damping. Our uniform sample of such intrinsically 
brightest pulses numbered 180 GPs with TOAs within 
the Fermi observing time. 

5. CORRELATION ANALYSIS 

The correlation analysis focused on two distinct tasks. 
The first one aimed to probe if the GP generation rate 
correlates with observed 7-ray ph otons. For the se cond, 
we investigated the hypothesis of iLvutikovl (120071 ) that 
the 7-ray photon flux increases during GPs. For both 
cases, we used simulated high-energy data sets with no 
assumed intrinsic correlation between the GPs and 7-ray 
photons to test the statistical level of correlation present 
in the real data. 

The simulations used the gtobssim software from the 
Fermi tools package. We used the latest version of in- 
strument response function, Pass6_v3 together with the 
same spacecraft /pulsar ephemeris as in real data analy- 
sis. 

We simulated the pulsar using the PulsarSpectrum li- 
brary, with the light curve , spectrum and int egral flux 
above 100 MeV taken from lAbdo et al.l (|2010D . The in- 
tegral flux was set to i^av = 2.09 x 10~^cm~^s~^ for the 
burst correlation analysis and varied from up to about 
a hundred i^av for the single-pulse correlation analysis 



(see I5.2|) . When we simulated zero flux from the Crab 
pulsar, we simply removed the pulsar from the list of 
simulated sources. 

We modeled the Crab Nebula as a point source (for 
the energy ranges in question its angular diameter is less 
than the Fermi r egion of inte r est), with the spectrum 
as determined in lAbdo et al.l ()2010[) and integral flux 
above 100 MeV of 9.8 x lO-'^cm-^s'^ For the Galactic 
and extragalactic backgrounds we used the "GalacticDif- 
fusc_v02" and "IsotropicDiffuse_v02" models. The simu- 
lated photon files were processed in the same way as the 
real data. 

5.1. Is GP rate correlated with single ^-ray photons? 

To test if 7-ray photons are correlated with the GP 
generation rate, we calculated the distribution of time 
lags between each photon and all GPs in that photon's 
GTI. The same procedure was applied to the simulated 
7-ray data sets, such that if there were any clustering of 
GPs around 7-ray photons (or with some time lag with 
respect to the 7-ray photons), it would be seen as a dis- 
crepancy between the real and simulated distributions 
of the high-energy data. Changing the bin size of the 
distribution makes it sensitive to different timescales of 
possible clustering of GPs. In this study, we tried a set 
of bin widths, starting from 10 seconds and increasing 
the width by 10 seconds up to 2 minutes. Two minutes 
corresponds to the smallest timescale of GP clumping 
(likely caused by interstellar scintillation), as shown in 
Fig. [5] On timescales less than 10 s, the Poisson noise 
due to a discrete number of time lag measurements be- 
comes too high. We performed 1000 simulation runs and 
contrasted the real-data distribution with the mean and 
95% percentile of all of the simulated data sets. 

Fig-Elshows the distribution of time lags between GPs 
and photons for all GPs and for a RISS-corrected sample 
of GPs for one particular bin width, namely, 2 min. In 
both cases the real data set lies all within the 95% per- 
centile of the simulations, indicating no apparent change 
in GP generation rate on 2-min timescales with any possi- 
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ble time lag within ±40 minutes (maximum GTI length) 
of the arrival of the 7-ray photons. All other bin widths, 
down to 10 s, gave the same result. 

5.2. Does ^ -ray flux change around single GPs? 

Another question of interest is whether the average 7- 
ray flux from the Crab pulsar increases during individual 
giant pulses, as predicted bv iLvutikovi (j2007i ). To inves- 
tigate that, we looked for the number of 7-ray photons 
in on-pulse emission windows around each CP. We per- 
formed a separate search for all GPs, looking for photons 
in a large window consisting of the main pulse, interpulse 
and bridge between them, and also for IP GPs only, limit- 
ing the correlation window to the interpulse phase range. 

If a photon was detected in a window around a GP, 
it was called a "match". For 10.5 hours of simultaneous 
observations, we detected only one such match: a photon 
with E = 403.7MeV was detected within 1.3ms of IP GP 
with peak flux density of 8.9 Jy. Nonetheless, knowing 
the observed number of matches iV = 1, the probability 
that the 7-ray flux during GPs is equal to some value Fq 
can be estimated with the simple Bayesian formula: 



p{F = Fo\N) 



pjFo) ■ p{N\F ^ Fq) 
Jo^'^^'' piFo) ■ piN\F = F^)dF^ 



(2) 



where p{Fq) is the prior distribution for Fq and p{N\F ~ 
Fq) is the likelihood, i.e. the probability to get the ob- 
served number of matches N if the pulsar 7-ray flux dur- 
ing GPs is equal to Fq. 

Since little is known about p{Fq), the prior distribution 
for Fq, we chose the prior to be uniform in a flux range 
from (the Crab pulsar turns off 7-ray emission during 
GPs) to 



F„ 



F^, 



observing timespan 
TVqp X length of window ' 



(3) 



where is the average pulsed 7-ray flux from the Crab 
pulsar, i^max corresponds to the hypothesis that all 7-ray 
photons from the Crab pulsar come during GPs. For our 
choice of windows, i^max ranged from 60 i^av (for on-pulse 
phase window) to 150 F^^ (for IP window). 

The likelihood p(A^ I F = Fq) was calculated by running 
simulations with different pulsed flux Fq and computing 
the fraction of runs with a number of matches N . The 
grid of trial flux values, in units of Fav, was as follows: 
from to 1 with the step of 0.25, from 1 to 20 with the 
step of 1 or 0.5, and then from 20 to 30 with the step of 
5. For both choices of correlation window the probability 
density went down to before 30 Fav Here we implicitly 
assumed that a higher flux outside selected windows does 
not influence the correlation within windows. 

Since the number of simulation runs for each trial Fq is 
finite, it leads to an uncertainty in estimating the likeli- 
hood. We estimated the statistical errors from the simu- 
lation using the following method. Suppose that for some 
value of Fq we have run n simulations with y successes 
(i.e. cases where the number of matches in the simula- 
tion equals the one obtained for real data, N). Then y /n 
defines the estimate of probability of success p, which is 
also the likelihood density p[N\F = Fq). More precisely, 
p\y has a Beta distribution, with mean [y+l)/ {n+1) and 

We adopted Up as an error 



of p due to limited numbers of simulations performed. 

However, there is another major source of uncer- 
tainty connected to the fact that we record discrete 
number of photons around the GPs. Since photon 
detection is very well described as a Poisson process 
(jRamanamurthv fc ThompsonI 119981 ). the error on de- 
tecting photons in a certain window around GPs will 
be a//V. In our case, for all windows we had A'' = 1, so to 
estimate the true value of the likelihood we should take 
into consideration also the likelihood curves for iV = 
and N = 2. These estimates for both windows are plot- 
ted in Fig. [21 left. The shaded region around each curve 
corresponds to iup, calculated by the above formula. 
For both correlation windows = 1 likelihood curves 
have maximum around Fg/Fav = 1, which means that, 
most probably, pulsed 7-ray flux does not change during 
GPs (no correlation) , or changes no more than few times 
(weak correlation). 

With our limited data set we cannot say anything more 
about the exact value of 7-ray flux during GPs, but we 
can place upper limits on it. On the grid of simulated 
fluxes Fi, one can convert the continuous formula for 
posterior probability density (eq. [2]) into a discrete one 
for the probability that pulsed flux around GPs is less 
than Fq: 



0.5-j:{p,+p,+i)iF,+i~F,) 



variance az. 



(y+l)(n-y+l) 
(n+3)(n+2) 



P{F<Fq\N)^ 

u.o ■ ; ,\Vi -V l[-ri+-\ — r.u 

(4) 

where pi = p{N\F = F^). 

To estimate errors in P{F < Fo|A^), we assumed that 
our uncertainty in pi due to a limited number of trials is 
much larger than the error from calculating the integral 
as a sum. As one can see in Fig. |6l left, this simplifica- 
tion is reasonable. Assuming all ap- are independent, the 
uncertainty in P{F < Fq\N) is determined by standard 
error propagation. 

In Fig. [HI we show the resulting probabilities that the 
7-ray flux from the Crab pulsar during GPs does not ex- 
ceed a given numb er of times the mean flux reported by 
lAbdo et al.l (j2010f) . Errors due to the limited number of 
simulation runs are plotted as errorbars, whereas those 
due to a discrete number of GP/photon matches are 
given by separate posterior probability curves for A^ = 
and N = 2. 

Obviously, the smaller the correlation window for a 
fixed observation timespan and the smaller the number 
of GPs in the sample, the larger Flnax, and the broader 
is the resulting posterior probability density. That is 
why for our data set we could obtain the posterior prob- 
ability densities only for correlation windows which in- 
cluded the IP, because most of GPs come within this 
phase range. For the main pulse GPs, p{F = Fo|A^) is 
very broad, having almost the same probability density 
up to ^ lOOFav. For the same reasons, the analysis on 
the sample of GPs corrected for refractive scintillation 
did not give any meaningful results. 

6. CONCLUSIONS 

No obvious correlation was found between Fermi pho- 
tons of energies > 100 MeV and radio giant pulses at 
the freguency of 8.9 GHz. No change in the Crab GP 
generation rate was found on timcscalcs from 10 to 120 s 
around 7-ray photons and with any possible lag within 
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FlG. 6. — Left: Likelihood, the probability of observing A'^ photons in a certain window around a GP if the 7-ray flux from the pulsar 
in this window were i^o. For each window size, we observed only one such match between GP and 7-ray photon, so Af = 1. Since photon 
arrival times are a Poisson process, the error on N is y/N = 1, thus likelihood curves for N = and N = 2 are also shown. The shaded 
region around each curve corresponds to uncertainty due to the limited number of simulation runs (see text for explanation). Both N = 1 
likelihood curves peak near Fq/Fhv = 1, so no or weak correlation between GPs and 7-ray photons is the most probable. Right: Posterior 
probability that 7-ray flux in a window around GP is less than Fq, given the observed number of matches A'^. Errors due to the limited 
number of simulation runs are plotted as errorbars, whereas those due to a discrete number of mat ches are given by s eparate posterior 
probability curves for N = and N = 2. i^av, the average pulsed 7-ray flux from Crab pulsar, is from lAbdo et aEI 12010 ). 



±40niin with respect to 7-ray photons. 

With 95% probabihty, the high energy flux of the Crab 
pulsar during GPs is less than 4 times the average 7- 
ray pulsed flux for the on-pulsc (MP-f-IP-t-bridgc between 
them) phase window. For IP GPs only, the 95% upper 
limit on 7-ray flux in the IP phase window is 12 times 
the average pulsed flux. If we consider the uncertainty 
due to discrete numbers of matches between photons and 
GPs, the 95% upper limits are 3—5.5 times the average 
pulsed flux for the pulsed emission window, and 8—16 for 
the IP window. 

A few explanations may be offered for the non- 
detection of GP-7-photon correlations. The most natu- 
ral is that production of GPs depends on non-stationary 
changing coherence conditions, which vary by a large de- 
gree even for similar magnetosphcric particle densities. 
Another possibility is that beaming in radio and at high 
energies are somewhat different, so that simultaneous 
GPs and 7-ray photons are emitted in different direc- 



tions. 

Overall, our results suggest that enhanced pair creation 
is not a dominant factor for GP occurrence, at least for 
high frequency IP GPs. However, our flux increase esti- 
mations are not on the leve l of a few percent, as in the 
work of iShearer et al.l ()2003[ ) at optical wavelengths. To 
reach that sensitivity we need more data, which will help 
push down the upper limit on flux during GPs and will 
make possible the analysis on subsamples of GPs, such 
as the brightest ones. Also, including radio frequencies 
below 4 GHz is potentially interesting not only for inves- 
tigating the correlation for MP GPs separately (MP GPs 
are much more common at lower frequencies), but also 
for re-doing the analysis for low-frequency IP GPs, since 
they might be generated by dif ferent physica l processes 
than the high-frequency IP GPs (;Mof fettlll997| ). AU these 
questions are being investigated with our ongoing radio 
observation campaign using the 42-ft telescope at the Jo- 
drell Bank Observatory (UK) and the 140-ft telescope at 
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the Green Bank Observatory (WV). 
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